A First Look at the Auriga-California Giant Molecular Cloud With Herschel and the CSO: Census of the Young Stellar Objects and the Dense Gas

We have mapped the Auriga/California molecular cloud with the Herschel PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech Submillimeter Observatory (CSO) with the eventual goal of quantifying the star formation and cloud structure in this Giant Molecular Cloud (GMC) that is comparable in size and mass to the Orion GMC, but which appears to be forming far fewer stars. We have tabulated 60 compact 70/160um sources that are likely pre-main-sequence objects and correlated those with Spitzer and WISE mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their properties. The most important result from this part of our study is that we find a modest number of additional compact young objects beyond those identified at shorter wavelengths with Spitzer. We also describe the dust column density and temperature structure derived from our photometric maps. The column density peaks at a few x 10^22 cm^-2 (N_H2) and is distributed in a clear filamentary structure along which nearly all the pre-main-sequence objects are found. We compare the YSO surface density to the gas column density and find a strong non-linear correlation between them. The dust temperature in the densest parts of the filaments drops to ~10K from values ~ 14--15K in the low density parts of the cloud. We also derive the cumulative mass fraction and probability density function of material in the cloud which we compare with similar data on other star-forming clouds.Comment: in press Astrophysical Journal, 201

Dissecting a hot molecular core: The case of G31.41+0.31

We made a detailed observational analysis of a well known hot molecular core lying in the high-mass star-forming region G31.41+0.31. This core is believed to contain deeply embedded massive stars and presents a velocity gradient that has been interpreted either as rotation or as expansion, depending on the authors. Our aim was to shed light on this question and possibly prepare the ground for higher resolution ALMA observations which could directly detect circumstellar disks around the embedded massive stars. Observations at sub-arcsecond resolution were performed with the Submillimeter Array in methyl cyanide, a typical hot molecular core tracer, and 12CO and 13CO, well known outflow tracers. We also obtained sensitive continuum maps at 1.3 mm. Our findings confirm the existence of a sharp velocity gradient across the core, but cannot confirm the existence of a bipolar outflow perpendicular to it. The improved angular resolution and sampling of the uv plane allow us to attain higher quality channel maps of the CH3CN lines with respect to previous studies and thus significantly improve our knowledge of the structure and kinematics of the hot molecular core. While no conclusive argument can rule out any of the two interpretations (rotation or expansion) proposed to explain the velocity gradient observed in the core, in our opinion the observational evidence collected so far indicates the rotating toroid as the most likely scenario. The outflow hypothesis appears less plausible, because the dynamical time scale is too short compared to that needed to form species such as CH3CN, and the mass loss and momentum rates estimated from our measurements appear too high.Comment: Astronomy and Astrophysics, in pres

Rotational Structure and Outflow in the Infrared Dark Cloud 18223-3

We examine an Infrared Dark Cloud (IRDC) at high spatial resolution as a means to study rotation, outflow, and infall at the onset of massive star formation. Submillimeter Array observations combined with IRAM 30 meter data in 12CO(2--1) reveal the outflow orientation in the IRDC 18223-3 region, and PdBI 3 mm observations confirm this orientation in other molecular species. The implication of the outflow's presence is that an accretion disk is feeding it, so using high density tracers such as C18O, N2H+, and CH3OH, we looked for indications of a velocity gradient perpendicular to the outflow direction. Surprisingly, this gradient turns out to be most apparent in CH3OH. The large size (28,000 AU) of the flattened rotating object detected indicates that this velocity gradient cannot be due solely to a disk, but rather from inward spiraling gas within which a Keplerian disk likely exists. From the outflow parameters, we derive properties of the source such as an outflow dynamical age of ~37,000 years, outflow mass of ~13 M_sun, and outflow energy of ~1.7 x 10^46 erg. While the outflow mass and energy are clearly consistent with a high-mass star forming region, the outflow dynamical age indicates a slightly more evolved evolutionary stage than previous spectral energy distribution (SED) modeling indicates. The calculated outflow properties reveal that this is truly a massive star in the making. We also present a model of the observed methanol velocity gradient. The rotational signatures can be modeled via rotationally infalling gas. These data present evidence for one of the youngest known outflow/infall/disk systems in massive star formation. A tentative evolutionary picture for massive disks is discussed.Comment: 11 pages, 9 figures. Accepted for publication in A&A. Figures 2,3,6, and 9 are available at higher resolution by email or in the journal publicatio

The Herschel view of the on-going star formation in the Vela-C molecular cloud

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.Comment: 9 pages, 7 figures, accepted by A&

A High Mass Dusty Disk Candidate: The Case of IRAS 18151-1208

Many questions remain regarding the properties of disks around massive prototstars. Here we present the observations of a high mass protostellar object including an elongated dust continuum structure perpendicular to the outflow. Submillimeter Array 230 GHz line and continuum observations of the high mass protostellar object IRAS 18151-1208 along with single dish IRAM 30m observations afford us high spatial resolution (0.8") as well as recovery of the extended emission that gets filtered out by the interferometer. The observations of 12CO confirm the outflow direction to be in the southeast-northwest direction, and the 1.3 mm continuum exhibits an elongation in the direction perpendicular to the outflow. We model the physical parameters of the elongated structure by simultaneously fitting the observed spectral energy distribution (SED) and the brightness profile along the major axis using the 3D Radiative Transfer code MC3D. Assuming a density profile similar to that of a low mass disk, we can also reproduce the observations of this high mass protostellar object. This is achieved by using the same density distribution and flaring parameters as were used in the low mass case, and scaling up the size parameters that successfully modeled the circumstellar disk of several T Tauri stars. We also calculate that a region within the inner 30 AU of such a high mass disk is stable under the Toomre criterion. While we do not rule out other scenarios, we show here that the observations in the high mass regime are consistent with a scaled up version of a low mass disk. Implications on high mass star formation are discussed.Comment: 14 pages, 11 figures, accepted for publication in Ap

A <i>Herschel</i> and BIMA study of the sequential star formation near the W 48A H II region

We present the results of Herschel HOBYS (Herschel imaging survey of OB Young Stellar objects) photometric mapping combined with Berkeley Illinois Maryland Association (BIMA) observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra-compact (UC) H II region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow and finally the NH2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UC H II region and discuss the Aquila supershell expansion as a major influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large-scale structures through kinematic information

Cluster-formation in the Rosette molecular cloud at the junctions of filaments

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.Comment: A&A Letter, in pres

Herschel Observations of a Potential Core Forming Clump: Perseus B1-E

We present continuum observations of the Perseus B1-E region from the Herschel Gould Belt Survey. These Herschel data reveal a loose grouping of substructures at 160 - 500 micron not seen in previous submillimetre observations. We measure temperature and column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the Green Bank Telescope. We find that the B1-E clump has a mass of ~ 100 solar masses and appears to be gravitationally bound. Furthermore, of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that B1-E may be forming a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our results suggest that B1-E may be influenced by a strong, localized magnetic field, but further observations are still required.Comment: 14 pages, 8 figures, published in A&A: Minor calibration correctio

Herschelobservations of the W3 GMC (II): clues to the formation of clusters of high-mass stars

The W3 giant molecular cloud is a prime target for investigating the formation of high-mass stars and clusters. This second study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further constrain the processes leading to the observed structures and stellar population. Probability density functions (PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as extinction AV. The shape of the PDF, typically represented with a lognormal function at low Av “breaking” to a power-law tail at high Av, is influenced by various processes including turbulence and selfgravity. The breaks can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Av(SF)» 6–10 mag) appears to shift to higher Av by stellar feedback, so that high-mass star-forming regions tend to have higher PDF breaks. A second break at Av> 50 mag traces structures formed or influenced by a dynamic process. Because such a process has been suggested to drive high-mass star formation in W3, this second break might then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency in a dense medium could result in a self-enhancing process, leading to more compression and favorable star formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be compatible with the “convergent constructive feedback” model introduced in our previous Herschel study